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Abstract 

We consider the stability of spatially homogeneous plane-wave space- 
times. We carry out a full analysis for plane-wave spacetimes in (4 + 1) 
dimensions, and find there are two cases to consider; what we call non- 
exceptional and exceptional. In the non-exceptional case the plane waves 
are stable to (spatially homogeneous) vacuum perturbations as well as a 
restricted set of matter perturbations. In the exceptional case we always 
find an instability. Also we consider the Milne universe in arbitrary di- 
mensions and find it is also stable provided the strong energy condition is 
satisfied. This implies that there exists an open set of stable plane-wave 
solutions in arbitrary dimensions. 

1 Introduction 

Finding solutions to Einstein's equations is a problem of great mathematical and 
physical interest. However, due to the fairly complicated nature of the Einstein 
equations, we usually impose symmetries in order to make the field equations 
more tractable. Some of the most successful schemes of symmetry reduction are 
the Bianchi models in (3+l)-dimensional cosmology 1,2,3,4 . For these models 
there exists a simply transitive group acting on the spatial hypersurfaces; hence, 
they are spatially homogeneous. For the spatially homogeneous models, the 
equations of motion reduce to a simple set of ordinary differential equations 
which makes it possible to study the evolution of these models in great detail. 
In particular, the asymptotic behaviour of all non-tilted Bianchi models with 
a perfect fluid have been determined [4,5,6 7 8 . Motivated by the success of 
the (3 + l)-dimcnsional Bianchi models, here, we shall investigate some of their 
higher-dimensional counterparts. 
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In (3 + l)-dimensional cosmology, algebraically special solutions usually play 
a particular role; they are equilibrium points, or sometimes even attractors, for 
more general classes of solutions. In particular the vacuum plane-wave solutions 
of dimension 3+1 are attractors for non-tilted non-inflationary models of Bianchi 
class B Hence, general solutions of Einstein's equations of Bianchi class 

B containing non-tilted non-inflationary perfect fluids can be approximated by 
vacuum plane-wave solutions at late times. For higher-dimensional models the 
situation is still unsettled as no stability analysis within this class of models has 
been done to date. 

In this paper we intend to take the first step towards this and analyse the 
(4+ l)-dimensional case. More precisely we look at spatially homogeneous plane 
waves, and find that they are stable to both perturbations in the geometry and 
also a restricted set of matter perturbations (which exclude inflationary fluids), 
provided a particular condition is satisfied; we call this the non-exceptional case. 
In the exceptional case, there exists an instability. The exceptional cases allows 
- due to the exact vanishing of one of the Einstein constraint equations - for 
an additional shear degree of freedom; it is this extra shear mode which turns 
out to be unstable. Also, as a precursor to future work, we show that the Milne 
universe is stable to perturbations in arbitrary dimensions provided the strong 
energy condition is satisfied. In particular this implies that in an arbitrary 
dimension, there exists an open set of stable plane-waves. 

Recently, of course, plane waves (and pp-waves) have featured heavily in the 
String Theory literature |9lll)llll| . Field theories in curved backgrounds usually 
assume that the back-reactions on the geometry can be neglected. However, 
the existence of an unstable mode would ultimately lead to the approximation 
being invalid. Therefore knowledge of the stability of such spacetimes in higher 
dimensions is crucial to such investigations, for studying quantum theories in 
such backgrounds inevitably leads to issues of back-reaction. We refer the reader 
to the works |12U13j where a different approach to analysing stability of plane- 
wave backgrounds, in a supergravity context, is taken. 

The outline of the paper is as follows. In section 2 we write down the 
relevant equations of motion for spatially homogeneous models. In section 3 we 
specialise to plane waves. Section 4 concerns the general stability analysis in 
4+1 dimensions. In section 5 we consider matter perturbations, and in section 
6 we analyse the stability of the (n + l)-dimensional Milne universe. Finally we 
provide a brief discussion. 

2 Equations of motion and constraints in n + 1 
dimensions for spatially homogeneous space- 
times 

In this section we recall some well known |3I14) results on spatially homogeneous 
spacetimes in n + 1 dimensions. Such spacetimes are defined as ones possessing 
a time-like vector n M hypersurface orthogonal, such that there exists a Lie group 
G acting simply transitively on the hypersurfaces. Note that this implies that 
the spacetime is of the form Rx £ with metric g = —dt 2 + h(t), where = idt)^ 
and h(t) is a left-invariant metric on S. For such spaces, n M is geodesic and the 
extrinsic curvature on the hypersurfaces reduces to = V^n,, = 9^ v . Now, 
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we can introduce an orthonormal frame on our spacetime, so 1 

g = vw"W = -dt 2 + 5 ab Lo a io\ (1) 

where we have picked w° = dt. The dual basis to will be denoted by {e^}, 
and we have 

[e^ej =7^„e p , (2) 

or equivalently dw M = T^ up uj u A lu p , where T pfJil/ = —j(-y pfU/ + j^p - Ivpn)- Note 
that orthogonality implies that 

7% = 0. (3) 
The remaining structure constants are 

7 a ob = -e a b + n\ (4) 

and 7 a bc = C a bc ; note O a f, is the local angular velocity of a set of Fermi- 
propagated axes with respect to {e a }, in the rest frame of an observer with 
velocity n M . The shear, a pu , is defined to be the traceless part of 9 pl/ , so 
@ij,v = Vi-iv + -9h^. 

The structure constants C a bc can be used to define the Lie algebra of the 
group G. The Bianchi models in 3 + 1 dimensions are constructed by considering 
all possible C a bc . Analogously to the 3+1 case we write these structure constants 
in terms of their traceless part D a bc , so D a ac = and 

C a bc = D a bc + 2S? b a c]l (5) 

where a b = -p^jC a ab . The Jacobi identity implies that D a bc a a — 0, as well as 
the dynamical equations, 

± ab + l dab _ (Q<* _ a f) ad = , (6) 
at n 

jD\ c + Ub\ c + 2(n d [b - af b )D\ ]d + (n a d - al)D d bd = 0. (7) 

So far the only assumption is spatial homogeneity. Now we impose the field 
equations. Einstein's equations, in a vacuum, are of course ^ n+1 ^Rfj,u = 0. 
These can be used in conjunction with standard results on hypersurfaces to give 
us our equations of motion. Firstly we have |14| . 



(n+1) 



R^rf = K 2 - K„ V K^ - V.KV^") + V M (n v V v n"), (8) 



which can be used to derive an (n + l)-dimensional version of Raychaudhuri's 
equation, namely, 



-6 2 + = 0, (9) 

n 



where '= n^W^. Another nice fact is |14|. 

^G^n" = I + K 2 - K^K^ , (10) 

iHere, Greek indices [fj,, u, ...) run over the full spacetime manifold, 0, l,...,n; lower case 
Latin indices (a, b, ...) run over the spatial hypersurfaces, 1, ...,n; and upper case Latin indices 
(A, B, ...) run over 1, n — 1. 
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which allows us to deduce the Friedmann constraint, 



i - - V = o^a^ - Mr. (ii) 



The < -™ +1 'i? a fc = h^h^ n+1 ^ components of the vacuum equations lead to the 
shear equations of motion, 2 



cr ab + 0a ab + {n >R {ab} = 0, (12) 

where A^ ab ^ of a tensor A ab means the trace-free part. One can calculate the 
curvature tensors on the hypersurface in terms of the structure constants. We 
get, 

{n) R ab = b ab - ^(D aM + D bad )a d - (n - l)a 2 h ab , (13) 
^R^b a a -n(n-l)a 2 , (14) 
b ab ee -±(2D c ad D cb d ~ D acd D b cd + 2D c ad D d bc ). (15) 

Finally the Gauss-Codazzi equation, 

D a K a b - D b K = (^i^nX (16) 

where D a is the metric connection of h, leads to d b C b ca — 0^C b ab = 0, which in 
turn leads to the constraint, 



na ac a a - a a b D b ac = 0. (17) 

Now wc will recast these equations in terms of expansion-normalised vari- 
ables following the (3 + l)-dimensional analysis @|. So we define H ab ee a ab /0, 
d b ac ee D b ac /9, n ab ee n ab /9, A b ee a b /6, ^K ab ee ^R ab /9 2 , and introduce a 
new time variable, r, by dr/dt ee 8. Then using Raychaudhuri's equation, the 
shear equations, the Friedmann constraint and the constraint implied by the 
Gauss-Codazzi equations, we find, 

Kb = (s 2 - ^) Sa h - 2^ Q fV - Wft (o6) , (18) 

1 - -) = S 2 - ( n) K, (19) 
n) 

n^ ac A a - ^d b ac = 0, (20) 

where ' = d/dr and E 2 = S o6 E o6 . 

At this point it is important to note that our equations of motion look 
like X' = F(X) where X is a column vector containing each of the independent 
dynamical variables. This is an autonomous system of differential equations, 

2 Note that & ab = -^a ab + 2a? a Q b -) c . 
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which lends itself to a straightforward stability analysis. Suppose X is an 
equilibrium point; i.e. 

F(X ) = 0. 

Consider perturbing about this, so X = Xo + SX. Then one finds, that we can 
rewrite the system as 

SX' = DF{X )SX + O(SX 2 ), 

where DF is the Jacobian matrix, i.e. the matrix of derivatives dFi/dXj. Hence, 
close to the equilibrium point, we can drop the quadratic terms in <5X so the 
stability depends only on the matrix Z)F(Xo). More precisely, the stability 
criterion for this system is that if the eigenvalues all have negative real parts, 
then we have stability; zero eigenvalues need more analysis and the existence of 
a sole eigenvalue with positive real part implies instability. 

In particular, if the equilibrium points are actually a space parameterised by 
m variables (as is the case in the subsequent analysis), then perturbing inside 
this space, i.e. F(X + SX) = 0, implies 

DF(X Q )SX = 0; 

namely, SX is an eigenvector of DF(X ) with zero eigenvalue, and thus, DF(X ) 
will possess at least m zero eigenvalues. 



3 Homogeneous plane-waves 

Plane-waves are solutions to Einstein's equations possessing a covariantly con- 
stant null vector. Now, define the null directions ry^ — -j= (u) n ± ui°), and 

similarly e± = ^-(e„ ± e ), so 

g = 2 V + r]- + S ab lo a uj b , (21) 

where A = 1,2, ...,n — 1. The existence of a null covariantly constant vector 
implies two things: firstly that it is a Killing vector and secondly that its metric 
dual (?(£,•) is closed. We can choose £ to be parallel to e_ so £ = /e_ where / 
is a function. Choosing a left-invariant frame (which we are free to do) implies 
that [£, eA\ — which gives, 

i a A0 =i a An, (22) 

e A (f) = 0. (23) 

Secondly, the metric dual to £, i.e. fr] + is closed so d(frj + ) = which gives, 

7° a6 + 7"a h -0, (24) 
C n ab = due to orthogonality, (25) 

together with a condition on / which we will not deal with. The wave front is 
required to be flat which tells us that there is an (n — l)-dimensional Abelian 
subalgebra and this forces, 

C A BC = 0. (26) 
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This implies that as — ^zjC a aB — ~^z\C n nB = 0, which means ajj = a&bn and 
upon using the evolution equation for a^, we get, 

a% = n\. (27) 

Thus we are left with only the following non-zero components, 



C A Bn = G A B , (28) 
1 A bo = C A Bn . (29) 

For every element 9 £ Mat (n — 1,R) we have a possible plane- wave solution. 
However there exists some redundancy; namely the group of rotations 0(n — 1) 
will simply redefine our orthonormal frame {w 4 }- Also, we are allowed an overall 
dilation, since the Einstein equations can be used to determine Q A a = (n — l)a. 
Thus the space of plane waves is isomorphic to Mat (n — 1, M.)/[0(n — 1) x R]; 
note the dimension of this space is \n(n — 1) — 1, see Appendix B. Now, the 
structure constants 7 a ob = — 9% + £l a b imply the following constraint, 

e A -n A B = e A B . (so) 

We will define Tab — 6(45), i.e the traceless part of 0. Thus we have, 

CAB = T( AB y (31) 
Finally, the expansion- normalised Tab will be denoted by Tab- 



4 General analysis in 4 + 1 dimensions 

In this section we will describe the main results of the paper. In Appendix 
A the equations of motion can be found for the 4+1 case. We provide the 
corresponding plane-wave solutions and find that these spacetimes are indeed 
stable to perturbations in the geometry in the non-exceptional case, whereas 
unstable in the exceptional case. We should make a remark on exactly what 
sort of perturbations we are dealing with: they correspond to perturbations in a 
subspace of spatially homogeneous models, namely ones for which the structure 
constants are C A B4 — Q A B with the rest vanishing 3 . Homogeneous plane- waves 
certainly belong to this class as we showed in the previous section, although it is 
not clear that all perturbations should lie in this space too. For now we assume 
so, however we will discuss this issue later. 

It should also be emphasised that when we refer to stability we mean in the 
technical sense; i.e. the perturbations are bounded, and not in the asymptotic 
sense; i.e. all perturbations tend to zero at late times. 

4.1 Parameterisation of independent variables 

First of all we note that the constraint l|20|l leads to ZU4 = in what we call the 
non-exceptional case (see Appendix A). We will use our freedom to rotate the 

3 A11 spatially homogeneous cosmological models of dimension 4+1 are classified in 1151 . 
We refer to this work for the other possibilities for the structure constants. 
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orthonormal frame to put Tab in upper triangular form; note this can always 
be done if it has real eigenvalues using an orthogonal matrix. The complex 
eigenvalue case can also be done, but this time one needs a unitary matrix (see 
Appendix B). We parameterise the variables by, 

S — 2Si S12 _ S13 

5 + Sx+V3S 2 £ 23 | (32) 

S13 £ 23 <5 + Si — y3S2 

-2n T12 T13 
r u; = I ti + %/3r 2 r 23 ) , (33) 
Ti - V3> 2 

where tab = Tab/ A. Note that 

Y 2 - 12S 2 + 6(£ 2 + E 2 ) + 2(£ 2 2 + £ 2 3 + £ 2 3 ). (34) 

The independent variables left after elimination of the constraints (see Appendix 

A) are Y*ab (for A ^ B), Si, E 2 and tab (for A < B), t x , t 2 ; the constraints 
are used to eliminate A and S. Thus we have 10 independent variables. Note 
that the definition of tab simplifies the equations of motion, in particular the 
diagonal components, which correspond to the group parameters (see Appendix 

B) , are constants of motion. Thus we see a further simplification and hence we 
are left with effectively only 8 independent variables for each group type. 

4.2 Exact solutions 

Here we write down the general plane- wave solution to the equations of motion. 
These are equilibrium points to the equations of motion. We find, 

n = P, t 2 = q, T12 = Qi, t 13 = Q 2 , t 23 = Q 3 , 

S12 = \AQi, E13 = \AQ 2 , S 23 = ^AQ 3 , 

where p = p 2 + q 2 + ^(Qf + Q 2 + Q 2 ), (35) 

and p, q, Q\, Q2, Q3 are constants. Their explicit metrics are given in |lfi| . how- 
ever, these are not needed in this analysis. Our family of plane waves is thus 
parameterised by five constants. Note that this means we expect the Jacobian 
matrix to have three zero eigenvalues, since we have already discarded the vari- 
ables Ti and t 2 on account of them being constants of motion. These three 
eigenvalues correspond to eigenvectors tangential to the set of exact plane-wave 
solutions. 

4.3 Stability 

The equations of motion form an autonomous system of first order differential 
equations, and as previously explained, one simply needs the eigenvalues of 
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the Jacobian matrix evaluated at the point we are perturbing about. Upon 
linearising our system we find the following eigenvalues: 



Ai =0, 

A 4 = - 

A 7 
A 8 



A 2 = 0, 
3 



A, 



2(p + 2)' 

3 ±i^Ql 



A 5 
A 6 



Ql 



Ql 



2(p + 2) 



Ql 



Ql 



(36) 



2(p + 2) 

Note that, ignoring the zero eigenvalues, we have 5 eigenvalues with mani- 
festly negative real parts (which, curiously are also all equal). Thus in the 
non-exceptional case we have stability. The zero eigenvalues were expected and 
correspond to moving in the space of plane-wave solutions. 



4.4 Exceptional case 

This occurs when T B A + Aa8 B is not invertible, and a ^ 0. Therefore the 
constraint leads to <J4b(t b a +4S b )A = which implies 2a^\A(2 — t\) = since 
t is upper triangular. Therefore we see the exceptional case corresponds to the 
particular group types where n = 2, which allows for an = f^i ^ in general. 
Let us restrict to a 2 4 = 034 = 0. The equation of motion for £41 is, 



£41 = (£ 2 -^-4S + 2£^ E41. (37) 

Note that the exact solution is unchanged, since S41 = together with the 
exact solution in the non-exceptional case solves all equations of motion. If we 
evaluate E 2 — | — 45 + 2Si at the exact plane wave solution, it is easy to show 
that it is always greater that Z/{p+ 2) 2 . Now, linearising the equation for £41 
we get, 



SE' 41 = AJS41, (38) 
3 

Therefore we have found an instability. 



where A > - ; (39) 



5 Matter perturbations 

Having shown that the homogeneous plane waves are stable to geometric per- 
turbations in (4 + 1) dimensions, it is desirable to test their sensitivity to the 
insertion of matter. It proves convenient to do so in (n + 1) dimensions. We 
begin with a stress-energy tensor for a perfect fluid, 

T flv = pn^ + Ph^. (40) 

From the Einstein equations, one finds 

R^nV = k + - —t-p) • (41) 

\n — 2 n — 1 / 
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We will consider matter satisfying the barotropic equation of state P = (7 — l)p. 
In this case the conservation of energy implies 



p + jp9 = 0. (42) 
Raychaudhuri's equation now becomes 

a ■ ^2,2 , ..J*!-* 



o» + Kpl^—f)=0. (43) 

We now define an expansion- normalised density, f2 = p- . We can then compute 
its derivative with respect to r, and using Ij43(l we obtain 

O' = -7 + - + 2E 2 + — . (44 

\ n n — 1 J 

To linear order, then, recalling that we are perturbing about the vacuum, we 
arrive at 

80! = SQ ( -7 + - + 2Eg ) . (45) 



Here £q corresponds to the vacuum plane-wave solution we are perturbing 
about. We immediately deduce that these spacetimes are stable against such 
matter perturbations provided that 

\ + < 7, (46) 

for n = 4. This follows from the fact that the linearised equation for 5£l only 
depends on 5tt and thus the eigenvalues of the Jacobian matrix are unaltered. 
Given that S 2 is a positive-definite quantity, it becomes clear demanding sta- 
bility in this model is slightly more restrictive than simply requiring the Strong 
Energy Condition, namely 7 > ^, to hold. Note that the case 7 = 0, which is 
equivalent to having a cosmological constant, is always unstable. On the other 
hand, if 7 = 1 as it does for dust, we can always find a stable plane wave. 



6 The Milne universe 

In this section we will use the equations set up in the first section to analyse 
the stability of the Milne universe in an arbitrary dimension. This spacetime 
can be defined as one with T A B = Y<ab = 0; note that this implies A 2 = 1/n 2 . 
Explicitly the metric, in horospherical coordinates, looks like, 

ds 2 = -dt 2 + t 2 [dz 2 + e 2z £W) 2 J • (47) 
If we perturb about this solution we get the linearised equations, 



71 1 _ 

6E' AB = 5E AB + (n - 1)A6T, AB) , (48) 

n 

Sf' A B = 0, (49) 
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where the ' = d/dr since we are working at linearised level. We see that this 
system of autonomous equations gives rise to a block diagonal Jacobian matrix 
whose eigenvalues are zero or negative. The zeros actually correspond to moving 
towards another plane-wave solution. Therefore the Milne universe is stable in 
the sense that the perturbations remain bounded. To see the origin of the zero 
eigenvalues, write, 

( tfZ) ) = ( o" o ) ( sf<Z) ) (50) 

where we have substituted A = 1/n which is true for the Milne universe (A = 
— 1/n can be excluded since we can pick its orientation). Therefore, the zero 
eigenvalue corresponds to 5£ ab = <^T(ab) > an d thus we see that we do not move 
off the surface T,ab = ^(AB) which contains our plane waves as we saw earlier. 

Interestingly, a simple argument allows us to say something about the stabil- 
ity of plane-waves in arbitrary dimensions from the above analysis. The Milne 
universe is the S 2 — » limit of the set of plane- wave solutions. Observing that 
the Jacobian matrix is a continuous function of the state space variables (which 
implies that its eigenvalues have to be so too) we can deduce that: provided that 
eq. \^.b\j holds, there will exist an open set of plane-wave solutions which are 
stable to spatially homogeneous perturbations. 

7 Discussion 

We have found that the stability analysis of vacuum spatially homogeneous 
plane-waves in 4+1 dimensions separate into two cases. For the non-exceptional 
case they are stable to (spatially homogeneous) perturbations in the geometry. 
They are also stable to matter perturbations provided 7 > 1/2 + 2Eq. This, of 
course, excludes inflationary fluids. 4 The exceptional case, however, is always 
unstable. This result is completely analogous to the (3 + l)-dimensional case 
where a similar instability for the plane- wave solution for the exceptional Bianchi 
type VF j, 9 [Q exists. 

It should be noted, however, if one allows for perturbations with respect to 
group types then two subtleties may occur. For some group types, even though 
they cannot be written as C A B4 = @ A B , there exist Lie algebra contractions for 
which the limit is an algebra of this type JH] . Also, for some special solutions an 
extra Killing vector may appear which may allow for simply transitive symmetry 
groups other than those considered here ^Hj- In any case, these possibilities only 
arise for a set of measure zero in the set of plane- wave solutions. 

Also, we found that the Milne universe in arbitrary dimensions is stable 
provided 7 > — , which is the strong energy condition for a perfect fluid. The 
implications of this are that in arbitrary dimensions, there always exists a neigh- 
bourhood of stable solutions (around the Milne solution). 

Of course, one can consider more general types of fluids than those consid- 
ered here; namely tilted (non-comoving) fluids, anisotropic stresses, dissipative 

4 The cosmological no-hair theorem tells us that in dimension 4+1 for 7 < 1/2 (e.g. infla- 
tionary type fluids) the late-time asymptotes are the flat FRW universe, which is consistent 
with our matter perturbation results. 
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fluids, p-form fluxes, etc. For example, if one takes a tilted perfect fluid, then 
generalising the arguments in |20j indicates that non-comoving velocities may 
have a significant effect on the late-time evolution of cosmological models. 5 An 
interesting extension of our work would be to include a systematic study of these 
more general fluids. 

Based on the above analysis, one is tempted to believe that all higher- 
dimensional plane-wave spacetimes, except for the exceptional models, within 
the class considered here, are stable for all spatially homogeneous vacuum per- 
turbations. We have already shown that some stable plane-wave spacetimes ex- 
ist in any dimension, and the stability has explicitly been shown for the (3+1)- 
and (4+l)-dimensional cases. However, whether this is true in all dimensions 
awaits to be seen. 
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A Equations of motion in 4 + 1 dimensions 

Here follows the differential equations and constraints used in the analysis of 
this paper. We will concern ourselves mainly with the expansion-normalised 
variables: 

v - a AB 

Tab = ?f- (51) 
S= 39 

Latin capital indices run from 1...3, the coordinates of the wave front. Further 
we have the following equations 

4(7 A4 a + (t%T B a = 

<7%{3a8 B + B A ) = (52) 
4cr44a - °b tB a = 0- 

In the non-exceptional case (defined as Q A B + 3a5g being invertible) we infer 
from the second line above a\ = = Recall that under our parameterisa- 
tion, (744 = —3s. The last of the above equations then reduces to, using l|51|) 

-12SA - T lAB f BA = 0. (53) 



5 A quick calculation shows that conservation of angular momentum in universes close to 
isotropy implies that the peculiar velocities are growing at late times if 7 > 1 + 1/n. 
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Under the useful field redefinition, tab = we can solve for S 

v ^ba 

l^ABT 



12 ' 

A 



(54) 



Although the function S appears in our parameterisation of the T, B , it will not 
appear on the right hand side on account of the tracelessness of t a b . 

The introduction of the r variable also simplifies the (first order) equations 
of motion, which we now produce for convenience: 



A! = (35 + £ 2 )A 
r/ = r 2 ' = 

rV = (-3£i - V3E 2 ) t\ + 2S 12 (^n + \/3r 2 ) 

t\' = - (3Sj - V3E 2 ) + 2S 13 (3n - V3r 2 ) - 2S 23 r 1 2 + 2£ 12 t 2 3 
r 2 3 ' = 2V3E 2 t 2 3 - 4\/3£ 2 3 t 2 

Sab' - (V - ^ S AB - X ab (55) 
where the symmetric matrix \ab is defined to be 

Xii = 6A 2 n - 2£ 12 2 - 2E 13 2 + \a 2 (t 12 2 + r 13 2 ) + ^ 2 (n 2 + r 2 2 ) 

X12 = ^— - — 3SiSi 2 — \/3S 2 Si 2 - 2Si 3 S 23 + -A 2 t\T\2 

. 2 ,2 ^13723 

H — r i2 r 2 + A — - — 

Xl3 = -^A 2 T 13 - 3SlSl 3 + V3Sl 3 i; 2 + ^A 2 TlTl 3 - ^A 2 Tl 3 T 2 

X22 = -3A 2 (n + V3r 2 ) + 2S 12 2 - 2£ 23 2 - %A 2 r 12 2 + -A 2 t 23 2 

o o 

+ ^ 2 (r 1 2 +r 2 2 ) + iA 2 r 13 2 
X23 - -^ 2r 23 + 2S 12 S 13 + 2V3S 23 S 2 - ^A 2 t 12 t 13 - V3A 2 t 23 t 2 
X33 = -3A 2 n + 3V3A 2 r 2 + 2S 13 2 + 2S 23 2 - ^A 2 (t 13 2 + r 23 2 ) 

+ \a 2 {t 2 + t 2 )+ 1 -A 2 t 12 2 . (56) 

In addition we have the Friedmann constraint, which allows us to eliminate A 
from the equations: 

H^GH- (57) 

where r 2 = t^V^j. This leads unambiguously to 

3 _ V 2 

^ - ¥zn- < 58 » 
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Note that the hrst equation of l|55|) was used to derive the subsequent equations 
for the t a b so it is included automatically. Thus, S and A have been removed 
and given t\ and T2 (the group type) we are left with eight dynamic variables. 

B Some results on n-dimensional Lie algebras 

Two Lie algebras, given by the structure constants C k ^ and C*^, are isomorphic 
(over R) if and only of there exist a matrix A € GL(n,M) such that 

C% = {^)\c l vq kP i X> j . (59) 
We are interested in the Lie algebras where 

c A BC = o, c A B „ = e A B . (60) 

For this class of algebras, the transformation l|59|) can be written in terms of 

© i B = ( A " 1 ) A A W' ( 61 ) 

in some suitable new orthogonal coordinates (marked with tildes), and <f> is 
a non-zero constant. Hence, the Lie algebras are determined in this case by 
equivalence classes of the map (in matrix form) 

0H^A _1 eA, A e GL(n- 1,R), <j> e R- {0}. (62) 

The 4> allows us to do a global rescaling of 0, while the conjugation leaves the 
characteristic equation, 

det(0 - Al) = 0, (63) 

invariant. Hence, the Lie algebras can be given in terms of the eigenvalues of 
modulo a constant rescaling. If \ are the eigenvalues, and |Ai| 7^ 0, say, then 
the (n — 2)-tuple, 

A = (ra-'w)' W£|A ' +i1 (64) 

will be invariant under the transformation (|62|1 . This means that different A 
yields non-isomorphic Lie algebras. If all eigenvalues are distinct, then A will 
uniquely determine the Lie algebra. However, in the degenerate cases, where 
two or more eigenvalues are equal, we have to use the Jordan canonical form 
to determine the Lie algebra. From this one sees that in general the space of 
distinct Lie algebras within the class given in ea. H6(J|l . is of dimension (n — 2). 

As an example, consider the 4-dimensional real indecomposiblc Lie alge- 
bras |21II22I23| . Here, the general types are the two two-parameter families A± g 5 
(three real eigenvalues) and A^ q 6 (one real and two complex conjugate eigenval- 
ues). The algebras A^i-Ai^ are the degenerate cases with non-trivial Jordan 
canonical form; A^i has three zero eigenvalues; A\ 2 has two equal and non-zero 
eigenvalues; A^z has two zero eigenvalues; and A^ has three equal non-zero 
eigenvalues. 

In our case, is a real matrix. This implies that the eigenvalues Aj are 
either real, or come in complex conjugate pairs. Due to the gauge symmetry 
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of determining the orientation of the frame {w" 4 }, we are allowed to perform a 
0{n — 1) rotation of the spatial frame to simplify the matrix 0. 

Let us first assume that all the eigenvalues A^ are real. Then for any © e 
Mat (n — 1,R) there exist an R 6 0(n — 1) such that R _1 0R is an upper 
triangular matrix. Hence, by a suitable choice of frame we can assume that 
is upper triangular with its eigenvalues along the diagonal: 



o a, ••• e 2 



(65) 



V • ■ ■ A„_i / 



Note that, as explained earlier, the ratio of the diagonal elements determine the 
group type. 

If some of the eigenvalues are complex, then it is clearly not possible to 
perform a real rotation so that the © is on upper triangular form. However, 
using a little trick, we can still put it into an upper triangular form. For any 
© 6 Mat (n — 1,R) there exist a U € U{n — 1) such that U^©U is an upper 
triangular matrix. Using unitary matrices we can put it into the form l|65[l with 
the eigenvalues along the diagonal. However, using a unitary "rotation" we have 
to be careful about interpreting the dynamical variables, u A B , T A and n A B . In 
matrix form, the unitary transformation implies the following transformation 
ona A B ,T A : 

itmUVU, T^U f TU, (66) 
and £l A B have to be reinterpreted as a u(n — l)-connection: 

n.->u t 4u + u t nu. (67) 

at 

Furthermore, the equations of motion and the constraints have to be amended 
similarly: 

a A B T B A ~ Tr(<rtT), a AB a AB ~ Tr(<r*tr), (68) 

etc. The equations of motion can now be solved in a similar manner. 

For the plane-wave spacetimes, we can use this to find the dimension of 
the space of plane-wave solutions. Basically, a plane-wave spacetime is given 
by specifying the upper-triangular trace-free matrix T. All other dynamical 
variables can be given, due to the Einstein equations and constraints, and the 
plane- wave criterion, in terms of T. Hence, there is an (n— 2)(n+l)/2-parameter 
family of homogeneous plane- wave solutions of which (n — 2) parameters specify 
the Lie algebra type. 
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